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ABSTRACT 

We derive the non-linear relation between the core flux F v of accretion powered jets at a given 
frequency and the mass M of the central compact object. For scale invariant jet models, the 
mathematical structure of the equations describing the synchrotron emission from jets enables 
us to cancel out the model dependent complications of jet dynamics, retaining only a simple, 
model independent algebraic relation between F v and M. This approach allows us to derive 
the F v —M relation for any accretion disk scenario that provides a set of input boundary con- 
ditions for the magnetic field and the relativistic particle pressure in the jet, such as standard 
and advection dominated accretion flow (ADAF) disk solutions. Surprisingly, the mass de- 
pendence of F v is very similar in different accretion scenarios. For typical fiat-spectrum core 
dominated radio jets and standard accretion scenarios we find F v ~ M 17 Z 12 . The 7-9 orders 
of magnitude difference in black hole mass between microquasars and AGN jets imply that 
AGN jets must be about 3-4- orders of magnitude more radio loud than microquasars, i.e., 
the ratio of radio to bolometric luminosity is much smaller in microquasars than in AGN jets. 
Because of the generality of these results, measurements of this F v — M dependence are a 
powerful probe of jet and accretion physics. We show how our analysis can be extended to de- 
rive a similar scaling relation between the accretion rate rh and F v for different accretion disk 
models. For radiatively inefficient accretion modes we find that the flat spectrum emission 
follows F v oc (M m) 17 ^ 12 . 

Key words: radiative mechanisms: non-thermal - galaxies: active - galaxies: nuclei - galax- 
ies: jets - x-rays: binaries - radio: continuum: general 



1 INTRODUCTION 

Relativistic jets are collimated outflows from the innermost regions 
of accretion disks around black holes and neutron stars. Not all 
accreting compact objects form jets, but when they do, the jet syn- 
chrotron radiation typically dominates the radio spectrum of the 
compact object. Such objects are called radio loud. 

Compact objects span 9 orders of magnitude in central mass: 
in active galactic nuclei (AGN), it ranges from M ~ 1O 6 M0 to 
M ~ few x 1O 9 M0, while Galactic X-ray binaries extend this 
range down to a few A/0. Yet, the jets formed by these objects ap- 
pear morphologically remarkably similar, and their core emission 
typically follows the same flat powerlaw spectrum. This suggests 
that the process of jet formation is universal, and that jets from 
supermassive black holes of different masses are not qualitatively 
different from each other and from jets in X-ray binaries, called 
microquasars. If this is so, we can compare jets from objects of 
different mass M (which is measurable dynamically) and accretion 
rate M (which is proportional to the accretion luminosity Lace) and 
determine how their observable characteristics change with M and 
M iSams et al J 199(1 . 

The most readily available observable parameter is the jet flux 



at a given frequency, F v . In this letter, we shall therefore derive, 
from theoretical arguments, the relationship between F v and M. 
As we will show, the mathematical structure of the expression for 
the jet synchrotron flux F v enables us to contract all the model 
dependent complications of jet physics into the formula for the ob- 
servable spectral index a and thus remove them from the relation 
between F v and M. For any observed value of a and for a set of 
boundary conditions delivered by accretion disk theory, we can thus 
formulate a model independent, non-linear relation between F v and 
M. 

A lot of effort has gone into searching for observational cor- 
relations between F„, M, and M. Such measurements are difficult 
because of numerous selection effects. Nevertheless, some obser- 
vational evidence of a non-l inear dependence between AGN ra- 
dio flux and black hole mass I Franceschini et al. 1998; Laor 2000; 
iMcLure &"b unlon 200ll: lLacv et alj |2001) exists in the recent lit- 
erature. However, other a uthors have found no such evidence 
I2002L I Woo & UrrvlEooll) . A systematic difference of radio loud- 
ness between neutron star and black hole X-ray binaries has also 
been suggested iFender & Kuulkersl200lh . 

Because of the large mass difference, any non-linearity be- 
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tween F v and M must be most apparent when comparing micro- 
quasars with AGN jets. And indeed, observations show that the ra- 
dio loudness parameter R — Lr/Luv, x-ray, defined as the ratio 
of radio luminosity (emitted by the jet) to UV/X-ray luminosity 
(emitted by the accretion disk) is much smaller for microquasars 
durin g outburst than it is for radio loud AGNs iFalcke & Biermannl 
Il99d) . In other words: the radio jet flux F v depends non-linearly 
on M. 

In the following sections, we will argue that this observational 
non-linearity is in not only consistent with, but required by the 
model independent F v -M relation that we will derive below. 



2 SCALE INVARIANT JETS 

The two fundamental parameters that determine the conditions in 
the inner accretion disk are the accretion rate M and the mass M 
of the central object. All length scales are proportional to the fun- 
damental scale of the compact object, r g oc M. The characteristic 
accretion rate of the disk is set by the Eddington rate MEdd, and 
for convenience we will define the dimensionless accretion rate 
rh = M/Mgad. All dynamically important variables are deter- 
mined by these two parameters. 

Because jets are formed in the inner disk, it is reasonable to 
assume that the conditions in the inner jet are set by the conditions 
in the inner disk, and thus, they will similarly depend on M and 
rri. Howeve r, it is possi ble that jets are powered by black hole spin 
extraction (Blandford & Znaiek 1977), in which case all jet vari- 
ables would also depend on the black hole spin parameter a. Thus, 
any dynamically important jet variable at the base of the jet will be 
determined by these three parameters. 

2.1 Dependence on M 

As mentioned in the introduction, observations suggest that the pro- 
cess of jet formation is universal, and that no qualitative difference 
exists between jets from objects of different mass. This morpholog- 
ical and spectral similarity of jets from objects with fundamentally 
different black hole masses suggests that jet formation and propa- 
gation might be a scale invariant processes 1 , i.e., that there is one 
relevant length scale in jet formation, which is r g , and that jet dy- 
namics are invariant under changes in this length scale 2 . 

1 The scale invariance assumed for the jet structure is only valid in the 
inner regions of the jet, where interactions with the environment are not im- 
portant. On large scales, where this interaction dominates the jet dynamics, 
additional parameters independent of the inner accretion disk enter, most 
prominently the external density and pressure. In this case, it is still pos- 
sible to write down extended scaling relations i Heinz 2002 ), but not in the 
form of eq. Q. However, since we restrict our analysis to the emission from 
the jet core, we can neglect these complications. 

2 While the influence of the spin parameter a on jet formation is not clear, 
it is important to note that a second length scale might be present in the 
process of jet formation, which is the light cylinder radius tq of the black 
hole spin. However, rn depends linearly on M and is thus a multiple of 
r g . The proportionality factor depends only on a, not on M, i.e., rn = 
/(a) r g . Thus, for a fixed spin parameter a, the only relevant length scale 
that changes upon variations in M, which we are primarily concerned with 
here, is r g . This is the length scale of relevance for changes in M, and we 
will henceforth assume that jet formation is invariant under changes in r g 
for otherwise fixed parameters, as suggested by the observational similarity 
of jets from very different M. If a does indeed have significant relevance 
for the process of jet formation, then variations in a will simply introduce 



Scale invariance implies that the spatial variation of important 
jet quantities (such as the shape of the jet, i.e., its lateral cross sec- 
tion, the orientation of magnetic field lines, the field strength, etc.) 
depends only on the dimensionless variable r/r g . Thus, a given 
variable / should be proportional to a function ip(r/r g ) which de- 
pends on r only through r/r g . In other words, we can scale a jet 
model for mass Mi by a length factor M2 /Mi and some spatially 
independent normalization factor and arrive at a jet model for mass 
Mi. 

In mathematical terms, this can be expressed as the condition 
that we can write any dynamically relevant quantity /, such as the 
magnetic field B(r), as the product of two decoupled functions: 

f(M,rh,a,r) = c/>/(M, rh, a) ijjf ^ — , rh, a^j 

= 4>f(M,m,a)ip f (x,m,a) (1) 

where r is the distance to the central engine measured along the 
jet, 4>J describes the dependence of / on the central engine mass 
M, and ipf describes the spatial dependence of / on the similarity 
variable % = r / r g f° r a given set of rh and a. Note that this is 
a requirement we place on the jet model, inspired by the observa- 
tional similarity between jets from different kinds of objects. Not 
all possible jet model must necessarily satisfy this relation. How- 
ever, those models that do satisfy it span an important sub-class of 
jet models and all of them will obey the rel ations derived below. 
One im portant example of such a model is the Blandfo rd & Konigll 
119791) model. 

The normalization functions c/> / reflect the dependence of the 
conditions at the base of the jet on the central mass M. Since jets 
are launched above accretion disks, it is natural to assume that these 
functions cA/ can be adopted from accretion disk models. 

For any geometric quantity, such as the jet diameter R(r), the 
direct proportionality of R to r g requires that <j)R tx M (where 
we have contracted the constant of proportionality into ipn for con- 
venience). As a dimensionless variable, it is reasonable to assume 
that the jet Lorentz factor F is entirely independent 3 of M. While 
measuring the bulk velocity of jets directly is impossible in most 
cases, the existing observational limits suggest that jets from mi- 
croquasars are not any more or less relativistic than their AGN 
counterparts, despite 9 orders of magnitude difference in M. We 
take this as sufficient evidence to assume in the following that 
r does not depend explicitly on mass, which allows us to write 
r(r) = Vt(X; a )> i- e -> 0r = 1. 

The state of the inner accretion disk depends on the accre- 
tion rate rh. In all standard accretion disk models, the fundamental 
quantities take on a rather simple scaling with the black hole mass 
M. For high accretion rates, where electron scattering becomes the 
dominant opacity source and radiation pressure dominates the ener- 
getics in the inner disk, the density and the pressure scale inversely 
with mass, q oc P ra d oc M" 1 fehakura & Sunvaevlll97rji . The 
magnetic field might be of the s ame order as the total pressure, 
and thus B oc M~ " 1/2 IShakura & Sunvae\ll973l) . Recent numeri- 
cal computations of the magneto-rotational instability in accretion 
disks seem to support this statement iBalbus & Hawlev 1998). The 
same scaling arises in advection dominated flows I Naravan & Yi 



intrinsic scatter to the relations derived below, as all relations will be derived 
for fixed, but arbitrary, a. 

3 Comparison with young stellar object jets indicates that the specific ve- 
locity Vv of jets is related to the orbital/escape speed of the innermost re- 
gions of the disk, which is independent of M for black holes. 
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1995), and in any scenario where dynamical terms dominate the 
cooling rate (e.g., convection or outflow dominated disks). Thus, 
we have <f)p = M~\ (j> e = M~\ and (f> B = M~ 1/2 . Observa- 
tions of X-ray binaries suggest that jet formation is linked to the 
so called low-hard state (Fender & Kuul kerlEoOll) . which is be- 
lieved to arise from an optically thin, geometrically thick accretion 
disk, which follows these scalings. We will therefore adopt (f>p — 
4>e = 4>b = M~ as our fiducial values. O nly if the innermost 
disk is of the standard gas pressure dominated Shakura & Sunvaevl 
Jl973l) type (which might be the case in low rh, low a accretion 
disks in AGNs. lFrank et alj2002h does this scaling differ slightly: 
here, P oc B 2 oc M" 9/1 ° and g oc M~ 7/w . 

Jets emit synchrotron radiation from a powerlaw distribution 
of electrons: 



dnjdrf — C7 p 



(2) 



within an energy range 7 m i n < 7 < 7 max . Here, 7 is the parti- 
cle Lorentz factor, C the normalisation constant, and p the power- 
law index. The production of powerlaw distributions is a universal 
property of diffusive shock acceleration, and we will assume that 
the fundamental powerlaw parameters p and 7 m i n are universal in 
relativistic jets as well. The observations from the optically thin 
part of jet spectra typically give p > 2. In the following, we will 
take p — 2 as our fiducial value for numerical examples. Since for 
spectra with p ^ 2 the high energy cutoff is dynamically unim- 
portant, we will not be concerned with its behavior. C is then di- 
rectly proportional to the pressure in relativistic particles, and we 
can once again write C = <f>c(M) t/'cKXi a )- ft is reasonable 
and customary to assume that the relativistic powerlaw particle dis- 
tribution is injected at some (unknown) fraction of equipartition 
with the magnetic field pressure, so C oc B 2 . Thus, for our fiducial 
values we have (j>c = M~ . 

The functions ip(x) can, in principle, take on rather compli- 
cated behavior, depending on the specific jet model. We will not 
be concerned with the detailed nature of ip, so long as they are 
mathematically well behaved (see i!3.2l for a definition of what this 
means). 



2.2 Dependence on rh 

While the main aim of this paper is to derive the scaling relation 
between jet radio flux and black hole mass, it is interesting to con- 
sider the dependence on other accretion disk parameters, namely rh 
(see also <i3.3i . The dependence of the fundamental disk parameters 
on m varies more significantly between different accretion models 
than the depende nce on M: 

jShakura & Sunvaevlll973l) showed that in radiation pressure 
supported disks, the total and magnetic pressure are independent 
of rh. Thus, the magnetic pressure B 2 and the relativistic particle 
pressure C in the jet are also independent of rh. The mass density 
in the disk, on the other hand, should follow g oc rh~ 2 M~ 1 , which 
might or might not affect the mass loading and thus the Lorentz 
factor of the jet. 

In mechanically cooled accretion disks (e.g., ADAFs), the 
pressure P and particle density g are directly proportional to rh. In 
jets launched from such flows, we thus have B 2 oc C oc mM _1 . 

We can derive the same scaling if we simply assume that the 
mechanical jet power W^t should be proportional to the disk lu- 
minosity Ldisk oc M — Mm: Since the jet power at injection is 
carried by internal energy, and since we can assume that the mag- 
netic field and the relativistic particle pressure are related to each 
other by some form of dissipation (e.g., reconnection, shocks), both 



C and B are related to the jet power by B 2 oc C oc W ]e t/R 2 c oc 
M/M 2 oc mM" 1 . As mentioned above, we chose this parametri- 
sation as our fiducial case. 

Finally, in standard gas-pressure dominated disks (however, 
still dominated by electron scattering, appropriate for the inner gas- 
pressure dominated disks of AGNs and X-ray binaries), the pres- 
sure follows P oc rh 4 / 5 M -9 / 10 , and thus, for the jet plasma we 
have B 2 oc C am' M~ 9/2 °, while the mass density follows 
pocm 2 / 5 M- 7 / 10 . 

After this excursion into the dependence of accretion disk and 
jet parameters on rh, we now proceed to investigate the radiative 
characteristics of self-similar jets. 



3 THE NON-LINEAR SCALING OF JET FLUX WITH 
BLACK HOLE MASS AND ACCRETION RATE 

3.1 Synchrotron emission from self similar jets 

The synchrotron self-absorption coefficient is 



A v C B 



(p+2)/2 -(p+4)/2 



(3) 



where A v is a p roportionality constant weakly dependent on p 
jRvbicki & Li ghtmarl 19791) . 

For ease of expression, we will present the following analysis 
in the case of a jet viewed from side on, however, extension to the 
case of arbitrary viewing angles is straight forward, and the result 
we derive is fully general. In the perpendicular case, the expression 
for t„ takes on a particularly simple form: 

R ict A p CB {p+2)/2 u~ [p+A)/2 



Rjct Ov 

A p M (f>c(M,rh,a 



p + 2 _ p + 4 

b(M, m,aj] 2 v 2 

p + 2 

■^R{x,m,a) 4) C (x,m,a) PU^.a)] 2 



= ®(M,m,a,v) $>(x,rh,a) 
where we define 

$ (M, rh, a, v) = M<f> c {M, rh, a) [4 



P + 2 

b(M, rh,a)\ 2 v 



(4) 



(5) 



®(x,m,a)=A p ipR(x,m,a)i>c(x,m,a)[4>B(x,™>,a)] 2 (6) 

The optically thin synchrotron emissivity for a powerlaw dis- 
tribution of electrons (well away from the lower and upper cutoff 
in the energy distribution) follows 

p+i p-i 
j u = J P CB 2 v 2 

= J p <j>c{M,m,a)[4> B {M,rh,a)] E ^ 1 



ipc(x,rn,a) [ipB(x,m,a)]~ 



(7) 



where J p is a con stant weakly dependent on p 
iRvbicki & Li ghtmarl 1 9791 . 

For simplicity, we will combine the dependence on the view- 
ing angle i9 due to Doppler beaming and optical depth effects into 
the function Because the viewing angle -3 and the Lorentz 

factor r are independent of M, it follows that must also 
be independent of M, which justifies this approach in what fol- 
lows. The jet surface brightness at a given frequency v is then 
S v ~ j v (1 — e~ T ")/a,,. The jet flux F v is then simply the 
surface integral over S„ : 



drR(r)S v {r) w C(#) / drR(r)j„{r) 



1 



a v {r) 
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ml dr [R(r)] 2 Mr) 1 ~ ' 



oc Ci^M^^c^v'^ 



r,,{r) 

dx4> R ^cij B 2 — 



p+i p-i 



oc M 0c S 2 v Q[$(M,rh,a,v),rn,a,-d\ 



(8) 



The integral O depends on M, and v only through the combination 
<3> from eq. J5j. 



3.2 The relation between F v and M 

From eq. JSJ, we can now work out the non-linear dependence 
of F v on the central engine mass M. The spectral index a = 
— din [Fv)j 91n (u) of the jet emission is given by 



91n(F„) 
d\n(v) 



p-1 9 In (9) 91n($) 



2 

p-1 



91n($) d\n{u) 
9 In (9) /p + 4^ 



(9a) 
(9b) 



2 91n($) 

Now taking the partial derivative of eq. with respect to M and 
substituting 9 ln(9)/91n(<I>) fromeq. I9bl . we can write 

p+i 



91n(F„) 
dln(M) 



= 3 + 



9 In 0c 9 In £ 



+ 



9 In (9) 91n($) 



2p- 



ain(M) 91n(M) 9 In ($) 9 In (M) 
13 + 2a 9 In {<j> B ) ( 2p + 3 + ap + 2a 



p + 4 
9 In (0 C ) 



+ 



91n(M) V p + 4 



91n(M) 
5 + 2a 



p- 



(10a) 



Quite generally, the functions 0c and 0b will be simple powers of 
71/ — for our fiducial assumptions, 0c = Af _1 and0s = M _1,/2 , 
and thus the index £m will simply be a constant: 



2p + 13 + 2a 
p + 4 

12 3 



2p + 3 + (p + 2)a 



0.33a 



5 + 2a 
p + 4 

(10b) 



where the approximate expressions assume p = 2. Thus, for any 
given set of rh, a, and #, F v will follow a simple powerlaw relation 
in M with powerlaw index £m 

F v OC ~ M l-42-0.33a 



(ID 



Variations in the other source parameters rh, a, the viscosity param- 
eter a v i sc , and ■& will only cause a mass independent scatter around 
this relation. 

Remarkably, this result is entirely independent of the functions 
ipf. Given a set of functions <f>f, which describe the dependence of 
the input conditions in the inner disk on M, and given an observed 
jet spectrum with spectral index a, eq. <1 It predicts the scaling of 
jet flux F v with M for any jet model that reproduces this spectral 
slope. The only assumptions that went into the derivation of this 
result are that a) the relevant parameters can be decomposed fol- 
lowing eq. Q, b) that the high (low) energy cutoffs in the spectrum 
are far above (below) the observed spectral band, and c) that the 
function 9 is analytic. This is what was meant when we required 
the functions ipf to be mathematically well behaved in ij2| 

Typically, the radio emission from core dominated jets follows 
a flat spectrum over many decades in frequency, i.e., a ~ 0. In this 
case, it follows for our fiducial parameters that the radio flux F r de- 
pends non-linearly on the mass to the £m = 17/12 ~ 1.42 power, 



once again, independ ent of the jet model, which manifests itself 
only through ipf. Falcke & Biermann 1 1995) based their adaptation 
of the original Blandford & Konial 1 1979) model on the assump- 
tion that B oc M~ 1//2 and C oc B 2 . And indeed, for their specific 
choice of ipf, they found £m = 17/12, which they already showed 
to be consistent with observations of flat spectrum radio jets from 
AGNs and microquasars (Falcke & Biermann 1996). 

As we mentioned before, the fiducial B 2 oc C oc M~ x scal- 
ing arises in a number of standard scenarios for the inner accretion 
disk: both in high efficiency, radiation pressure dominated inner 
disks and in low efficiency ADAFs. The value of £ = 17/12 — a/3 
is therefore a very general result, which depends only weakly on 
the spectral index a. In jets that are launched from standard gas 
pressure dominated disks that extend all the way to the innermost 
stable orbit, the change in the scaling to fl = tpc = M _9//1 ° 
leads to a change in the Mass index: £m = [143 + 22p — a(14 + 
9p)]/[20(p + 4)] ~ 1.56 — 0.23a, which is even more non-linear 
than the standard value of £m ~ 1.42 — 0.33a from eq. jlObt . If 
the magnetic field responsible for spin extraction from black holes 
is supported by or anchored in the inner disk, the same considera- 
tions might hold for Blandford-Znajek (1977) jets. 

It is worth noting that this analysis holds even for the case of 
jets composed of discrete ejections or internal shocks, if we define 
F v as the time averaged flux or the peak flux. In fact, because the 
derivation of eqs. 181 1 Oat did not assume any specific jet-like geom- 
etry, they hold for any synchrotron emitting plasma with powerlaw 
spectrum if the source parameters can be described by eq. Q. 



3.3 The relation between F v and rh 

An interesting feature of the derivation of eq. ilOai is that it is mod- 
ular: any fundamental accretion disk parameter that enters into the 
dynamical description of the jet in an invariant fashion following 
eq. such that it only appears in the functions 0/ and not in ipf 
leads to such a relation. For example, if we can separate the depen- 
dence of any dynamical quantity / on the accretion rate rh into a 
function 0/ so that f(M, rh, a, r) — <p(M,rh, a) ip(x,a), we can 
derive a non-linear relation between F v and rh of the form 4 



dln(F v ) _ 9 In (0 S ) / 2p + 3 + a(p + 2) 



9 In (rh) 9 In (rh) V p + 4 

+ 91n(0 c ) f 5 + 2a 



a In (m) V p + 4 ' 

following the same derivation as in eq. 

For our fiducial assumption 0c oc <p 2 B oc rh (from ADAF type 
accretion, or the Ansatz Wj c t oc Ldisk) we get 

where the approximate expressions assume p = 2. Note that for 
flat spectrum sources (a = 0), the dependence on rh is the same as 
that on M, as found by Sa ms et alJ Il996f) . 

Using the accretion disk scaling relations discussed in i|2.2l 



4 Note that the assumption that the jet Lorentz factor Tjct is indepen- 
dent of rh, which is implicit in combining line-of-sight effects into an rh- 
independent function is not necessarily given. In the case where a 

strong dependence of rj ct on m arises, the complications introduced by 
Doppler beaming will introduce an rh-dependent scatter in the F v -m rela- 
tion, which could skew the distribution away from the mean scaling index 
£m expected from eq. I12al . 
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Table 1. The dependence of B and C on Af and m and the scaling indices 
5m an d 5m for different accretion modes (rows 1-3) and for the Ansatz 
that the mechanical jet luminosity Wj ct should be proportional to the disk 
power Ldisk ( row 4), assuming p = 2. 



we have presented the powerlaw indices for the jet scaling relations 
with mass M and accretion rate rh expected for these different ac- 
cretion modes in Table Q The last two columns in this table show 
the scaling indices £m and £ m , such that 

F„ocM^ M rh^ (13) 

If jet production by black hole spin extraction is invariant un- 
der changes in a in the sense of eq. Q, i.e., if it depends only 
trivially on a, we can write a relation between F v and a by simply 
replacing rh by a in eq. J12at 

3.4 Optically thin vs. optically thick emission 

Synchrotron self-absorption is stronger at lower frequencies. Thus, 
at high frequencies the jet must be optically thin even at the loca- 
tion where it is injected. Since eqs. J 1 Oat and 1 1 2al were derived 
without any restrictions on t v , we can use them to infer the scal- 
ing of radiation at high, optically thin frequencies as well, as long 
as the other assumptions made above still hold. The only assump- 
tion which might be violated is that radiative cooling of the elec- 
tron spectrum is negligible. Since modifications by radiative cool- 
ing should be visible as spectral breaks or cutoffs, we will continue 
to neglect them here, assuming that a spectral band can be chosen 
where the spectrum is optically thin yet unaffected by cooling. 

If the jet is injected at a distance Xi from the black hole (ex- 
pressed in dimensionless units), then the frequency v T= \ at which 
the jet spectrum becomes optically thin is given by eq. J4) by de- 
manding that T v (xi, M, rh,a) = 1: 



M<t>c4>T^(Xi, 



(m4>c4> 



I- I 2 
2 

B 



(14) 



Above v T , the spectrum is optically thin and has a spectral index of 
a = (as can be seen from eq. {8j by setting t u <C 1). Here 
and below, we will denote optically thin values by a tilde. E.g., Fo 
is the flux at an optically thin frequency v. 

To derive the scaling indices £m and £ m in the optically thin 
case, we could go through the same arguments as in eqs. jSJ through 
I12bi . now imposing that = <3> 1 I' <C 1. More easily, however, 
we can derive £m and £m by simply replacing a in eqs. HOai and 
I12ai by the optically thin value a = (p — l)/2. This gives 



6 



3 + 



d In (0c) , p+ldln(0 fl ) 



dln(M) 



+ 



2 91n(M) 



(15) 



where the expression in square brackets is valid for 0c oc B oc 
M" 1 and p = 2. 

For optically thin jets withp ~ 2. ISams et alj il9 961 suggest 
that the observed brightness temperature Tb i0 b s in microquasar and 
AGN jets decreases with bolometric luminosity as ( Mm) -0 ' 76 . 
The optically thin radio flux then goes as Fp l0 bs oc iifetTb.obe oc 



M Tb. obs M ' . Thus, the observations give ^M.obs = 1-24, 
which coincides nicely with the theoretical value of £m = 1.25. 
For the scaling of optically thin flux Ft with m, we find 



d In (<t>c 
d In (rh) 



+ 



.+ 19 In (0 B ) 
2 ain(m) 



(16) 



% oc rh 



where the expression in square brackets is valid for (j>c oc 
and p = 2. 

For fixed M, a change in rh results in a change to both the 
optically thick flux F v and the optically thin flux Fa. We can thus 
relate these changes to see how the optically thick flux varies as a 
function of the optically thin flux: 



d\n{Fu 



d In {Fa) 



(17a) 



where is given in eq. 1 1 2at and £ m in eq. 1161 . 

If we impose a unique relation between 0c and 0s (the most 
reasonable assumption here is a fixed fraction of equipartition be- 
tween relativistic particles and B-field at the base of the jet, which 
reproduces our fiducial assumption of 0c oc </>%), we can actually 
rewrite £p as 



p + 4 



5 + 2a + 



1 - p + 2a 

, , o | jj 4>c 

+ P + Z d In <t> B 



(17b) 



where rh is now only an implicit parameter. This implies that the 
relation holds for any variation in the jet parameters (whether it is 
caused by a change in rh or any other parameter) as long as it does 
not affect the geometry or dimensions of the jet (since both the 
functions tpf and Rj ct were kept fixed when deriving eq. I17bl ). 
For the fiducial assumption of 0c oc S , eq. I17bl reduces to 



4 + 2a + 



2a 



p + 5 



17 + 8a 
21 



(17c) 



where the expression in square brackets holds forp = 2. 

Thus, we have a relationship between optically thick and op- 
tically thin flux under the condition of fixed M of the form 



F 



(17+8t»)/21 



(18) 



which is remarkably close to the observed correlation between the 
flat spectrum radi o flux and the X-ray flux o bserved in the Galactic 
source GX 339-4 JCorbel et alJ200fll2003ll of ^ obs ~ 0.7. 

Once again, all of these relations arise independently of the 
specific jet model, so long as it produces a sp ectral index of a in the 
optically thick part of the spectrum. In fact, (Markoff et alJ J2003I) 
find exactly the sam e result for £ g, when applying their jet model 
iMarkoff etall200lh to GX 339-4 as eq. fTTcl for the case of a = 
and p = 2. 



4 DISCUSSION 

4.1 Observational Consequences 

Because the relations of eq. 1 1 Oat and eq. <12at are model- 
independent (they only depend on the boundary conditions at the 
base of the jet, not on ipf), measurements of £m, £m, and £p cannot 
be used to distinguish between different jet models. However, the 
generality of this result makes such measurements an even stronger 
probe of the underlying nature of jet physics: 
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• Observational confirmation of eqs. HOai and I12ai would 
prove that jet formation is scale invariant. On the other hand, if ob- 
servations can rule out a any correlation which can be described by 
these equations, this would argue strongly against scale invariance. 

• Measuring the values of f; ai and £™ would provide diagnostics 
of the conditions in the inner disk and at the base of the jet, i.e., 
measuring £m and £m could be used to put limits on cj>f, 

• Because in low efficiency accretion the accretion rate cannot 
be measured directly, it might not be possible to establish a direct 
observational correlation between F v and rh. However, if the above 
relations hold, any correlation of the optically thick, flat spectrum 
jet radio emission with emission at higher frequencies, could be 
used to constrain the high energy emission processes (e.g., optically 
thin jet emission or bremsstrahlung from an ADAF). 

• Measuring the residual spread of F v around the predicted re- 
lation could provide a handle on the relative importance of orienta- 
tion effects, and thus to measure the mean jet Lorentz factor Y. 



4.2 Conclusions 

We have derived the non-linear relation between the observed jet 
flux at a given frequency F v and the black hole mass M. For scale 
invariant jets, the nature of the expression for the jet synchrotron 
emission makes it possible to contract all the model dependence 
into the observable spectral index a. Thus, for any observed value 
of a, the derived F v -M relation is now model independent — any 
jet model that produces the observed jet spectrum automatically 
satisfies this relation. Given a prescription of the input conditions 
at the base of the jet, provided by accretion disk theory, we can thus 
predict the scaling of jet flux with M, Most accretion scenarios pro- 
duce a scaling relation of the form F v <x M 17 / 12- "/ 3 . Thus, for 
the optically thick flat spectrum radio emission from core dominate 
jets we find F v ~ Af 17//12 , while for optically thin emission with 
a ~ 0.5, we find F v ~ M 5 ^ 4 . Due to the large range in black hole 
mass, this non-linearity makes AGN jets much more radio loud than 
microquasar jets. 

This analysis can be extended to any fundamental accretion 
parameter (e.g., accretion rate or black hole spin) if jet dynamics are 
invariant with respect to changes in this parameter. For example, for 
ADAF-like boundary conditions at the base of the jet, the scaling 
with accretion rate m follows F v ~ (M m) 17 ^ 12+2a ^ 3 , and in the 
flat spectrum case of a = 0, the dependence on rh and M is the 
same: F v oc M 17 ' 12 . Because this result is model independent, 
observational measurements of the non-linear scaling of F v with 
M and rh are powerful probes of the behavior of the underlying 
accretion flows and of the nature of the energy and matter supply 
to jets from compact objects. 

Clearly, the physics of jet formation is extremely complicated. 
For example, we still do not understand the nature of the radio- 
loudness dichotomy in AGNs, and even though they contain black 
holes of similar mass, GRS 1915+105 is much more active in the 
radio band than Cyg X-l. Nevertheless, independent of the com- 
plicated physics of jet formation, the arguments presented in this 
paper show that the radio loudness of jets increases with increas- 
ing black hole mass., and thus, that the radio emission from mi- 
croquasars should be a much smaller fraction of their bolometric 
luminosity than that of radio loud AGNs. 

We would like to thank Eugene Churazov, Tiziana DiMattero, 
Torsten Ensslin, Heino Falcke, Rob Fender, Sera Markoff, and An- 
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